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ABSTRACT 

We  report  the  discovery  of  a  nearby  old  halo  white  dwarf  (WD)  candidate  from  the  Sloan  Digital  Sky 
Survey  (SDSS).  SDSS  J1 10217.48+41 1315.4  has  a  proper  motion  of  1"75  yr-1  and  redder  optical  colors 
than  all  other  known  featureless  (type  DC)  WDs.  We  present  SDSS  imaging  and  spectroscopy  of  this  object, 
along  with  near-infrared  photometry  obtained  at  the  United  Kingdom  Infrared  Telescope  (UKIRT).  Fitting 
its  photometry  with  up-to-date  model  atmospheres,  we  find  that  its  overall  spectral  energy  distribution  is 
fit  reasonably  well  with  a  pure-hydrogen  composition  and  Teff  ~  3800  K  (assuming  logg  =  8).  This 
temperature  and  gravity  would  place  this  WD  at  35  pc  from  the  Sun  with  a  tangential  velocity  of  290  km  s-1 
and  space  velocities  consistent  with  halo  membership;  furthermore,  its  combined  main-sequence  and  WD 
cooling  age  would  be  ^11  Gyr.  However,  if  this  object  is  a  massive  WD,  it  could  be  a  younger  object 
with  a  thick  disk  origin.  Whatever  its  origin,  the  optical  colors  of  this  object  are  redder  than  predicted  by 
any  current  pure-hydrogen,  pure-helium,  or  mixed  hydrogen-helium  atmospheric  model,  indicating  that  there 
remain  problems  in  our  understanding  of  the  complicated  physics  of  the  dense  atmospheres  of  cool  WDs. 

Key  words:  stars:  individual  (SDSS  J1 10217.48+41 1315.4)  -  white  dwarfs 
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1.  INTRODUCTION 

There  is  considerable  interest  in  white  dwarf  (WD)  stellar 
remnants  in  the  halo  of  our  galaxy.  Since  halo  stars  are  gen¬ 
erally  older  than  disk  stars,  the  oldest  halo  WDs  should  be 
older  and  cooler  than  the  oldest  disk  WDs.  Such  objects  are 
of  interest  for  studying  the  age  distribution  of  halo  stars  and 
for  testing  our  understanding  of  stellar  atmospheres,  since  their 
spectra  can  depart  significantly  from  simple  blackbodies  (see, 
e.g.,  Hansen  1998;  Saumon  &  Jacobson  1999;  Kowalski  2006a). 
Halo  WD  candidates  may  be  identifiable  as  faint  objects  with 
high  proper  motions  and  unusual  colors  (e.g.,  Ducourant  et  al. 
2007).  Some  halo  WD  candidates  have  indeed  been  found,  but 
most  of  them  have  relatively  warm  temperatures,  indicating  that 
they  are  relatively  young  WDs  (e.g.,  Bergeron  et  al.  2005 ;  Lepine 
et  al.  2005).  To  date,  the  coolest  known  probable  halo  WD  is  WD 
0346+246,  with  T  ~  3800  K  (Bergeron  2001).  The  coolest  halo 
WD  candidate ,  whose  velocities  are  consistent  with  either  a  halo 
or  a  thick  disk  origin,  is  SDSS  J122048.65+091412.1  (Gates 
et  al.  2004).  The  latter  object  is  one  of  a  handful  of  ul¬ 
tracool  WDs  (reff  <  3800  K)  whose  optical  spectra  show 
collision-induced  absorption  (CIA)  from  H2.  CIA  causes  WDs 
to  exhibit  increasingly  bluer  colors  at  increasingly  shorter 
wavelengths  at  Te ff  <  5000  K  (see,  e.g.,  Bergeron  et  al. 
2005). 

Here  we  report  the  identification  of  a  candidate  old  halo  WD 
with  red  optical  colors  whose  spectrum  was  obtained  by  the 
Sloan  Digital  Sky  Survey  (SDSS;  York  et  al.  2000).  The  SDSS 
used  a  drift- scanning  imaging  camera  (Gunn  et  al.  1998)  on  a 
2.5  m  telescope  (Gunn  et  al.  2006)  to  image  ~  104  deg2  of  sky  on 


the  SDSS  ugriz  magnitude  system  (Fukugita  et  al.  1996;  Hogg 
et  al.  2001 ;  Smith  et  al.  2002;  Pier  et  al.  2003;  Ivezic  et  al.  2004; 
Tucker  et  al.  2006).  Two  multi- fiber,  double  spectrographs  are 
being  used  to  obtain  R  ~  2100  spectra  for  ~106  galaxies  and 
~105  quasar  candidates  (Stoughton  et  al.  2002).  As  discussed 
in  Richards  et  al.  (2002),  most  quasar  candidates  are  targeted 
for  spectroscopy  because  they  are  outliers  from  the  stellar  locus. 
Spectroscopy  of  such  targets  provides  data  not  just  on  quasars, 
but  also  on  objects  with  colors  different  from  those  of  the  stellar 
locus,  such  as  the  unusual  WD  presented  herein. 

2.  IMAGES  AND  PHOTOMETRY 

SDSS  J1 10217.48+41 13 15.4  (hereafter  J1102+4113)9  was 
noticed  during  visual  inspection  of  all  SDSS  spectra  from 
Data  Release  Six  (Adelman-McCarthy  et  al.  2008)  classified 
as  unknown  by  the  SDSS  pipeline.  It  stood  out  in  the  SDSS 
Catalog  Archive  Server  as  having  a  high  proper  motion  of 
fia  =  —105.0  ±3.5  mas  yr_1  and  /x5  =  —1750 ±3.5  mas  yr_1, 
computed  as  described  in  Munn  et  al.  (2004)  by  combining 
astrometry  from  the  SDSS  and  from  the  USNO-Bl.O  catalog 
(Monet  et  al.  2003).  J 1 102+41 13  is  present  in  seven  Palo- 
mar  Observatory  Sky  Survey  (POSS)  plates  and  is  present 
in  one  published  SDSS  observation.  A  selection  of  these  im¬ 
ages  is  shown  in  Figure  1.  J1 102+41 13  is  cataloged  in  USNO- 
Bl.O  as  1312-0217226,  with  proper  motions  of  pa  =  — 106  ± 


9  Due  to  its  high  proper  motion,  the  SDSS  coordinates  for  this  object  are  in 
Equinox  J2000  and  Epoch  2003.3.  The  United  States  Naval  Observatory 
USNO-B  1.0  coordinates  for  this  object  in  Equinox  J2000  and  Epoch  2000  are 
R.A.  =  11:02:17.50  and  decl.  =  +41:13:21.52. 
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I  POSS-II  ,R  50948, _ ,  ‘  ~  r  •  r  SDSS  ugriz  52754 

Figure  1.  Imaging  of  J1 102+41 13  (just  below  each  arrow)  at  six  different  epochs,  centered  on  its  position  in  the  epoch  of  the  SDSS  image  at  lower  right.  The  image 
source,  equivalent  passband,  and  Modified  Julian  Date  (MJD)  are  indicated  in  each  image.  J1 102+41 13  moves  from  the  top  middle  to  the  center  of  the  images  over 
time.  Each  image  is  5'  x  5'  in  size,  with  north  up  and  east  at  left. 

(A  color  version  of  this  figure  is  available  in  the  online  journal) 


2  mas  yr-1  and  /jl$  =  —1744  d=  1  mas  yr-1  based  on  five  pho¬ 
tographic  epochs  only.  We  adopt  a  total  proper  motion  of  /x  = 
1.75  arcsec  yr-1 .  In  Galactic  coordinates,  J 1 102+411 3  is  located 
above  the  Galactic  plane  near  anticenter  (/,  b  =  174°, 63. 5°). 
Its  proper  motion  is  predominantly  in  the  +/  direction.  Taken 
together,  these  findings  mean  that  its  velocity  parallel  to  the 
Galactic  plane  is  less  than  the  Sun’s  (see  Section  4.1). 

Despite  its  presence  in  the  USNO-Bl.O  catalog,  this  object 
has  not  been  reported  in  the  literature  as  a  high-proper-motion 
star.  No  sources  are  listed  in  the  SIMBAD  or  NASA/IPAC  Ex- 
tragalactic  Database  (NED)  databases  within  V  of  its  position. 
J1 102+41 13  is  bright  enough  that  it  could  have  been  found  by 
Luyten  (1974).  It  was  missed  during  examination  of  USNO- 
B  high-proper-motion  stars  by  Levine  (2005)  because  it  is 
0.3  mag  fainter  than  the  limit  of  R  <  18.0  used  in  that  study. 


It  should  also  have  been  found  in  the  LSPM-N  catalog  (Lepine 
et  al.  2003;  Lepine  &  Shara  2005),  although  it  is  near  the  lower 
magnitude  limit  and  upper  proper-motion  limit  of  that  catalog.  It 
may  have  been  missed  by  previous  searches  because  the  image 
on  the  POSS-I E  plate  appears  double  (Figure  1,  top  left  panel). 
This  is  likely  due  to  a  plate  flaw,  since  the  second  component  of 
the  putative  double  is  narrower  than  the  point- spread  function 
(PSF). 

Available  photometry  from  published  observations  of 
J 1 102+41 13  is  given  in  Table  1.  The  object  is  consistently  faint 
in  blue  passbands.  There  is  a  faint  J- band  feature  in  the  Two  Mi¬ 
cron  All  Sky  Survey  (2MASS;  Skrutskie  et  al.  2006)  within  2"  of 
the  expected  position  of  J1 102+41 13  at  that  epoch.  To  measure 
the  flux  in  this  feature,  we  retrieved  the  2MASS  Atlas  images 
covering  this  object,  measured  3"  radius  aperture  magnitudes 


00 


Table  1 

Photometry  of  SDSS  J1 10217.48+41 1315.4 


Source 

MJD 

U  ±  Gu 

r  ±  <jr 

i  ±  Gi 

Z  ±  crz 

J±oj 

H±gh 

K±gk 

u- g 

8-z 

POSS1 

35183.7938 

20.59 

POSS1 

35183.8321 

18.58 

POSS2 

47540.9117 

18.67 

POSS2 

48001.6992 

20.26 

POSS2 

48294.8804 

Unavailable 

POSS2 

50094.9138 

18.41 

2MASS 

50912.8346 

17.78  ±  0.76 

17.37  ±0.92 

17.52  ±  1.81 

POSS2 

50948.6813 

Unavailable 

SDSS 

52754.1402 

23.01  ±  0.48 

20.20  ±  0.02 

18.76  ±0.02 

18.21  ±0.02 

17.93  ±  0.02 

2.81  ±0.48 

2.27  ±  0.03 

UKIRT 

54456.1192 

17.24  ±  0.02 

17.33  ±  0.02 

17.34  ±  0.03 

Notes.  Optical  magnitudes  are  on  the  AB  system  (but  see  text),  and  2MASS  and  UKIRT  magnitudes  on  the  Vega  system.  2MASS  magnitudes  use  the  2MASS  filter  system  (Skrutskie  et  al.  2006)  and  UKIRT 
magnitudes  use  the  MKO  filter  system  (Ghinassi  et  al.  2002).  POSS  magnitudes  have  been  recalibrated  to  the  indicated  SDSS  filter  (Sesar  et  al.  2006).  SDSS  magnitudes  are  PSF  magnitudes,  but  have  been 
converted  from  asinh  magnitudes  (Lupton  et  al.  1999)  to  traditional  magnitudes,  using  the  information  in  Table  21  of  Stoughton  et  al.  (2002).  For  this  object,  the  difference  is  significant  only  in  the  w-band. 
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Figure  2.  Full  SDSS  spectrum  of  J1 102+41 13,  smoothed  by  a  7  pixel  boxcar, 
plotted  as  F\  (in  units  of  10~17  erg  s-1  cm-2  A-1)  vs.  wavelength  in  A.  The 
uncertainty  at  each  pixel  is  plotted  along  the  bottom. 

in  /,  H ,  and  Ks  at  the  position  of  the  potential  /-band  de¬ 
tection,  and  calculated  appropriate  magnitude  uncertainties.10 
These  measurements  are  included  in  Table  1,  but  their  large 
uncertainties  make  them  of  limited  use.  To  obtain  better  near- 
infrared  (near-IR)  constraints,  we  obtained  JHK  photometry  of 
J1 102+41 13  on  2007  December  21  using  the  United  Kingdom 
Infrared  Telescope  (UKIRT)  Fast-Track  Imager  (UFTI;  Roche 
et  al.  2003)  in  a  service  mode.  The  data  were  reduced  in  the 
standard  fashion  and  the  photometry  was  calibrated  using  ob¬ 
servations  of  UKIRT  faint  standard  #130  (Leggett  et  al.  2006). 
These  measurements  are  also  included  in  Table  1 . 

3.  SPECTRA 

Based  on  initial  (“TARGET”)  photometric  reductions  of 
SDSS  imaging,  J1 102+41 13  was  targeted  as  a  high-redshift 
quasar  candidate  (target  flag  QSO_HIZ).  A  single  SDSS  spec¬ 
trum  was  obtained  on  Modified  Julian  Date  (MJD)  53046  on 
the  SDSS  spectroscopic  plate  no.  1437  and  fiber  no.  428.  In 
Figure  2,  we  present  the  full  SDSS  spectrum  of  J1 102+41 13.  Its 
spectrum  is  red  and  featureless  (in  particular,  there  is  no  sign  of 
B aimer  absorption)  except  for  a  possible  broad  emission  feature 
at  5750-5900  A.  This  feature  is  near  one  end  of  the  wavelength 
region  5800-6150  A  where  both  the  blue  and  the  red  SDSS 
spectrographs  record  flux  from  an  object.  To  investigate  the  re¬ 
ality  of  this  feature,  we  examined  the  four  individual  blue  and 
red  exposures  which  were  all  combined  to  produce  the  final 
weighted  average  SDSS  spectrum  shown  in  Figure  2.  The  broad 
feature  is  present  in  each  of  the  four  individual  blue  exposures, 
but  not  in  any  of  the  four  individual  red  exposures.  In  fact,  the 
red  exposures  all  show  a  dip  in  flux  below  6050  A.  In  Figure  3 
we  present  the  average  spectra  of  this  object  obtained  with  the 
blue  and  red  SDSS  spectrographs,  shown  in  blue  and  red  respec¬ 
tively,  along  with  the  =bla  uncertainties  for  both  spectra,  shown 


5600  5800  6000  6200  6400 

Wavelength  in  Angstroms 

Figure  3.  Average  spectra  of  J1 102+41 13  from  the  blue-wavelength  SDSS 
spectrograph  (in  blue)  and  the  red-wavelength  SDSS  spectrograph  (in  red), 
along  with  ±la  error  bars  (in  gray).  The  discrepant  flux  levels  in  the  two 
spectra  between  5800  and  6000  A  casts  doubt  on  the  reality  of  the  emission 
feature  at  5750-5900  A  in  the  combined  spectrum. 

(A  color  version  of  this  figure  is  available  in  the  online  journal) 

in  gray.  The  flux  levels  at  5800-6000  A  in  the  two  spectra  do  not 
agree  within  the  uncertainties.  This  disagreement  is  indicative 
of  some  problem  with  the  SDSS  spectrum  in  this  wavelength 
range.  We  conclude  that  the  broad  feature  at  5750-5900  A  is 
an  artifact,11 * * V  although  an  independent  spectrum  would  still  be 
worth  obtaining  to  verify  that  conclusion. 

There  is  no  other  statistically  significant  absorption  or  emis¬ 
sion  feature  in  the  spectrum.  In  particular,  there  is  no  sign  of  Ha 
absorption  to  a  3 o  limit  of  ~0.8  A.  (The  smoothed  spectrum 
shows  a  dip  in  flux  at  6522  A,  but  it  is  due  to  two  noise  spikes 
which  are  narrower  than  the  instrumental  resolution.) 

4.  ANALYSIS  AND  COMPARISON  WITH  MODELS 

We  classify  J 1 102+41 13  as  a  DC  WD,  since  its  optical 
spectrum  has  no  robust  features.  The  featureless  spectrum  of 
J 1 102+41 13  is  similar  to  other  cool  DC  WDs.  At  optical  wave¬ 
lengths,  the  reddest  of  these  are  WD  0346+246  (Oppenheimer 
et  al.  2001b),  GD  392B  (Farihi  2004),  WD  1247+550  (Liebert 
et  al.  1988),  and  WD  1310-472  (Bergeron  et  al.  2001).  These 
four  stars  all  have  B  —  V  ~  1.4  and  1.39  <  V  —  I  <  1.46. 12 
Only  one  of  these  four  WDs  has  a  measured  g  —  i  color,  but 
the  small  dispersion  in  their  V  —  I  colors  means  that  there  will 
be  a  small  dispersion  in  their  g  —  i  colors.  Thus,  comparing  the 
g  —  i  =  1.67  measured  for  WD  1247+550  with  the  g  —  i  = 
1.99  measured  for  J 1 102+41 13,  we  conclude  that  J 1 102+41 1 3 
is  ~0.3  mag  redder  in  g  —  i  than  the  reddest  previously 


10  The  calculation  of  the  magnitude  uncertainties  followed  the  procedure  in 
section  VI.8.a.ii  of  the  2MASS  Explanatory  Supplement  at 
http://www.ipac.caltech.edu/2mass/releases/allsky/doc/explsup.html. 


11  Some  pixels  at  nearby  wavelengths  are  flagged  as  potentially  being 
untrustworthy,  with  some  or  all  of  the  flags  NEARBADPIXEL,  LOWFLAT, 

S C ATTEREDLIGHT,  B ADFLUXF ACTOR,  and  BADSKYCHI  (Stoughton 
et  al.  2002).  Empirically,  we  have  found  that  the  grow  radius  around  such 
untrustworthy  pixels  is,  on  rare  occasions,  not  large  enough  to  flag  all 
apparently  problematic  pixels  in  SDSS  spectra.  This  spectrum  may  be  one 
such  case. 

12  Two  WDs  reported  to  have  V  -  I  >  1.5  by  Oswalt  et  al.  (1996)  in  fact  have 

V  —  I  <  1.0  (Bergeron  et  al.  2001;  H.  C.  Harris,  unpublished). 
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known  DC  WDs.  The  only  other  non-magnetic  WDs13  known 
to  have  optical  colors  approaching  those  of  J1 102+4113  are  DZ 
WDs  with  extremely  strong  calcium  and  sodium  absorption: 
WD  2251-070  (Liebert  et  al.  1988),  WD  J2356-209  (Oppen- 
heimer  et  al.  2001a),  and  SDSS  J133001. 13+643523.8  (Harris  et 
al.  2003).  In  contrast,  J1 102+41 13  shows  no  evidence  of  absorp¬ 
tion  from  Na  i  AA5891,  5897,  Ca  i  A.4227,  or  Ca  n  AJU8500, 
8544,  8664.  A  cool  atmosphere  with  an  extremely  low  metal 
abundance  seems  the  best  explanation  for  its  red  color. 

For  a  surface  temperature  as  low  as  that  indicated  by  the 
red  color  of  J 1 102+41 13,  both  pure-hydrogen  and  pure-helium 
model  atmospheres  are  expected  to  be  essentially  featureless 
in  the  optical.  In  cool,  pure-H  atmospheres,  most  H  is  in  the 
form  of  H2  and  most  atomic  H  is  in  the  ground  state,  leading 
to  negligible  B aimer  absorption  (which  would  be  extremely 
pressure  broadened  in  any  case).  In  cool,  pure-He  atmospheres, 
most  He  is  neutral  and  the  lower  levels  of  optical  transitions  of 
He  i  are  not  populated  for  Te ff  <  12,000  K,  leading  to  highly 
optically  transparent  atmospheres.  Nonetheless,  there  can  be 
detectable  differences  between  the  spectra  of  hydrogen-  and 
helium-dominated  cool  WDs,  particularly  in  the  near-IR  where 
CIA  is  affected  by  the  composition. 

To  constrain  the  atmospheric  composition  of  J 1 102+41 13, 
we  fit  model  WD  atmospheres  to  the  UKIRT  photometry  and 
to  the  SDSS  photometry  converted  to  an  AB  system.  The 
SDSS  magnitudes  are  already  on  an  AB  system  to  within  ~1% 
in  g  and  r,  while  in  other  bands  we  applied  the  corrections 
wab  =  wsdss  —  0.04,  z‘ab  =  z'sdss+0-02,  andzAB  =  zsdss+0.03 
(Abazajian  et  al.  2004;  Eisenstein  et  al.  2006;  see  also  Holberg 
2007;  J.  A.  Holtzman  et  al.  2008,  in  preparation).  To  account 
for  ~1%  uncertainties  in  these  corrections,  we  increased  the 
uncertainties  from  2%  to  3%  for  each  SDSS  magnitude  besides 
u  in  the  fitting. 

In  Figure  4  we  plot  the  SDSS  spectrum  of  J 1 102+41 13,  the 
SDSS  and  UKIRT  photometry  and  synthetic  photometry  from 
fits  discussed  below,  all  in  units  of  flux  density  per  unit  frequency 
versus  wavelength.  For  illustrative  purposes,  the  spectrum  has 
been  scaled  upwards  by  0^6  to  match  the  optical  photometry. 
Such  scaling  is  needed  because  not  all  the  light  from  the  object 
is  captured  by  the  spectroscopic  fiber.  The  scaling  is  greater 
than  the  typical  SDSS  value  of  0^35  (Adelman-McCarthy  et  al. 
2008),  probably  because  the  fiber  placement  did  not  account  for 
the  l".  4  proper  motion  of  J1 102+4113  between  the  imaging  and 
spectroscopic  epochs. 

4.1.  Comparison  with  Pure-Hydrogen  Models 

We  first  fit  pure-hydrogen  atmosphere  models  (Kowalski 
2006b;  Kowalski  2007;  Kowalski  &  Saumon  2006)  to  the 
photometry  of  J 1 102+41 13  (excluding  the  u  band,  whose 
uncertainties  are  too  large  to  be  useful),  assuming  a  typical  WD 
gravity  of  log  g  =  8  (e.g.,  Fontaine  et  al.  2001).  The  resulting 
fit  yields  Te ff  =  3830  K.  Synthetic  photometry  from  this  model 
is  shown  as  the  blue  triangles  in  Figure  4.  The  fit  is  reasonable 
given  the  model  uncertainties,  although  the  formal  x2  is  quite 
poor:  x2  =  111  for  four  degrees  of  freedom  (v  =  n  —  1  =  6 
minus  Te ff  minus  a  normalization  factor),  with  about  ~75%  of 
the  x2  signal  coming  from  the  discrepancy  at  g.  The  red  optical 
colors  are  produced  in  the  far  red  wing  of  Lya  (perturbed  at 
high  density),  while  the  dips  in  the  near-IR  spectrum  are  due  to 
H2-H2  CIA. 

13  No  known  magnetic  WD  has  a  red,  featureless  spectrum  like  that  of 
J1 102+41 13.  Nonetheless,  polarimetric  observations  would  be  useful  to  ensure 
that  J 1102+41 13  is  not  some  sort  of  unusual  magnetic  WD. 


Figure  4.  Observed  spectrum  (plotted  as  Fv  in  units  of  erg  s-1  cm-2  Hz-1  vs. 
wavelength  in  microns),  observed  photometry,  and  synthetic  model  photometry 
of  J1 102+41 13.  The  SDSS  spectrum  is  shown  as  the  dark  gray  line  and  the 
associated  uncertainties  as  the  light  gray  line.  The  emission  feature  we  believe  to 
be  spurious  has  been  interpolated  over.  The  observed  photometry  of  J 1 1 02+4 113 
is  shown  as  the  black  error  bars.  Blue  triangles  are  synthetic  photometry  from 
the  best- fit  pure-H  model.  Red  squares  are  synthetic  photometry  from  a  pure-He 
model  fit  to  the  optical  data  only. 

(A  color  version  of  this  figure  is  available  in  the  online  journal) 

While  J 1 102+41 13  is  (F22  redder  in  g  —  z  than  any  pure-H 
model,  it  is  an  open  question  how  well  such  models  reproduce 
the  relevant  physics  at  Te ff  ~  3800  K.  Hydrogen  in  the 
atmosphere  of  a  Tc ff  ~  3800  K  WD  approaches  a  density 
where  nearby  molecules  are  strongly  correlated  and  where  the 
refractive  index  is  significantly  greater  than  unity.  Furthermore, 
there  are  still  uncertainties  about  the  reliability  of  H2  CIA 
opacities  in  WDs.  Observationally,  there  are  few  WDs  known 
near  that  temperature  with  which  comparisons  with  models  can 
be  made.  At  Te&  <  3800  K,  CIA  will  affect  the  optical  colors 
of  WDs  and  their  g  —  z  colors  will  become  bluer.  However, 
given  all  the  physical  effects  discussed  above,  which  may  not  be 
handled  correctly  in  the  models,  it  is  not  certain  how  optically 
red  pure-H  atmospheres  become  before  they  turn  bluer.  In  short, 
the  models  could  be  sufficiently  in  error  to  explain  the  formally 
poor  fit  of  the  pure-H  model  predictions  to  the  J 1 102+4113  data. 

The  distance  to  the  WD  in  our  pure-H  model  fit  is  35  pc, 
making  the  tangential  velocity  utan  =  290  km  s_1.  The  cooling 
time  to  reach  Te ff  =  3830  K  is  ^9.6  Gyr,  assuming  a  typical  WD 
mass  of  0.6  MQ  with  a  core  of  equal  parts  C  and  O  (Fontaine  et  al. 
2001).  Including  ~1  Gyr  for  its  progenitor’s  lifetime  (Fontaine 
et  al.  2001),  for  a  pure-H  model  J 1 102+41 13  is  the  remnant 
of  a  star  that  formed  ^11  Gyr  ago.  Alternatively,  at  the  same 
temperature,  if  log  g  =  9  (9.5),  then  J1 102+41 13  would  have 
the  following  approximate  properties:  M  =  1.11  (1.36)  M0, 
d  =  17  (10)  pc,  utan  =  140  (83)  km  s-1,  cooling  time  7  (3)  Gyr, 
and  progenitor  lifetime  1  Gyr.  If  J 1 1 02+4 113  has  log  g  =  7 . 5 , 
then  it  would  be  at  d  =  42  pc  with  M  =  0.35 M0  and  its  total 
age  would  be  greater  than  the  age  of  the  universe  unless  it  is  an 
unresolved  double  degenerate  or  a  product  of  common-envelope 
binary  star  evolution  (Fontaine  et  al.  2001). 
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Figure  5.  Synthetic  spectra  of  pure-H,  pure-He,  and  mixed  H/He  models,  all 
assuming  log  g  =  8.  Fits  to  all  data  yield  Teff  =  3830  K  for  pure-H  (solid  green) 
and  Teff  =  5381  K  for  pure-He  (solid  blue).  Fits  to  optical  data  only  yield  Te\j  = 
3450  K  for  pure-H  (dotted  green)  and  Teff  =  3360  K  for  pure-He  (dotted  blue). 
Helium-dominated  models  with  Te\j  =  3500  K  are  shown  by  the  red  lines,  with 
hydrogen  contents  as  follows:  log  H/He  =  —3  (solid),  —4  (dash-dotted),  —5 
(dotted),  and  -5.5  (dashed). 

We  can  calculate  the  components  of  the  object’s  space 
velocity  relative  to  the  Sun,  U,  V,  and  W,  which  are  positive  in 
the  directions  of  Galactic  center,  Galactic  rotation,  and  the  north 
Galactic  pole,  respectively.  If  we  assume  zero  radial  velocity, 
we  find  ( U,V,W )  =  (63,-280,46)  km  s-1  for  d  =  35  pc, 
(30,-140,22)  km  s"1  for  d  =  17  pc,  and  (18,-80,13)  km  s"1 
for  d  =  10  pc.  Chiba  &  Beers  (2000)  give  (([/>,  (V),  (W))  = 
(17  =b  141,  —187  zb  106,  —5  zb  94)  km  s-1  for  the  halo  and 
(4  zb  46,  -20  zb  50,  -3  zb  35)  km  s”1  for  the  thick  disk.  Thus, 
the  space  velocities  of  J1 102+41 13  are  most  consistent  with 
those  of  the  halo.  However,  if  J1 102+41 13  is  at  d  <  13.5  pc, 
it  could  be  a  young,  thick  disk  WD  with  an  unusually  (but  not 
unprecedentedly)  high  mass.  A  parallax  measurement  is  needed 
to  discriminate  between  these  possibilities  and  to  pin  down  the 
mass  of  J1 102+41 13. 

4.2.  Comparison  with  Pure -Helium  Models 

We  next  fit  pure-He  models  (Kowalski  &  Saumon  2004, 2006) 
to  the  data  (again  assuming  log  g  =  8),  but  were  not  able  to  find 
a  realistic  fit.  If  the  near-IR  photometry  is  omitted,  a  good  fit 
is  found  with  Te ff  =  3360  K.  Synthetic  photometry  from  this 
model  is  shown  as  the  red  squares  in  Figure  4;  the  model  near- 
IR  fluxes  greatly  overpredict  the  observed  fluxes.  Thus,  a  good  fit 
cannot  be  found  to  the  combined  optical  and  near-IR  photometry 
of  J 1 102+41 13  with  a  pure-He  model:  if  the  temperature  is 
increased  to  move  the  peak  of  the  fit  to  shorter  wavelengths, 
the  slope  of  the  optical  fit  will  deviate  unacceptably  from  the 
observations  (see  the  blue  lines  in  Figure  5). 

4.3.  Comparison  with  Mixed  Hydrogen-Helium  Models 

We  also  consider  mixed  H/He  models,  as  even  a  tiny  amount 
of  hydrogen  in  a  helium-dominated  atmosphere  (H/He  > 
10-10)  can  induce  sufficient  H2— He  CIA  to  cause  a  significant 
near-IR  flux  deficit  (see,  e.g.,  Figure  5  of  Bergeron  &  Leggett 
2002).  The  results  of  normalizing  mixed  H/He  models  with 
Te ff  =  3500  K  to  the  photometry  are  shown  in  Figure  5;  for 
slightly  higher  or  lower  Te ff  (zb500  K),  the  overall  results  are 
similar.  None  of  the  mixed  H/He  models  yields  a  better  match  to 
the  data  than  the  pure-hydrogen  model,  although  a  model  with 
H/He  =  10-5  fits  at  wavelengths  A.  <  1.5  /xm. 


In  summary,  the  poor  fit  of  pure-He  models  to  the  photometry 
means  that  hydrogen  is  certainly  there  in  the  atmosphere  of 
this  WD.  However,  the  H/He  ratio  cannot  be  pinned  down 
with  current  models,  due  to  uncertainties  in  how  accurately  the 
models  treat  the  complicated  physics  involved  (refraction,  non¬ 
ideal  equation  of  state,  non-ideal  chemistry,  perturbed  Lyman 
series  absorption,  etc.;  see  Kowalski  2006a).  Both  H/He  ^  1 
and  H/He  <Cl  are  possible. 

4.4.  Modeled  and  Observed  Fluxes  and  Colors 

If  only  observational  uncertainties  are  considered,  no  current 
WD  model  atmosphere  provides  a  formally  acceptable  fit  to  the 
photometry  of  J1 102+41 13.  To  guide  future  models,  it  is  worth 
examining  at  what  wavelengths  the  discrepancies  arise. 

In  both  Figures  4  and  5,  it  can  be  seen  that  J 1 102+41 13 
is  redder  at  optical  wavelengths  than  most  model  fits.  To 
even  approximately  match  both  the  red  optical  colors  and 
the  depression  of  the  near-IR  flux  (relative  to  a  blackbody 
of  the  same  temperature)  requires  some  hydrogen  in  its 
atmosphere. 

In  terms  of  u  —  g  versus  g  —  z,  J 1 102+41 13  matches  neither 
pure-H  models  nor  pure-He  models  (see  Table  1  and  Figure  4(b) 
of  Kowalski  &  Saumon  2006).  It  is  redder  in  g  —  z  than  pure-H 
models  at  any  u  —  g,  but  is  bluer  in  u  —  g  than  the  pure-He 
models  at  its  g  —  z  color.  Because  of  the  large  uncertainty  in  its 
u  —  g  color,  J 1 102+4 113  is  only  1.7a  redder  in  u  —  g  than  the 
pure-He  sequence  of  Kowalski  &  Saumon  (2006). 

5.  SUMMARY 

SDSS  J 1 10217.48+41 13 15.4  is  a  cool  WD  with  either  a  pure- 
hydrogen  atmosphere  at  Te  ff  ~  3830Kora  mixed  hydrogen- 
helium  atmosphere  with  H/He  ~  10-5  and  reff  ~  3500  K. 
In  either  case,  its  distance  is  expected  to  be  d  <  40  pc.  At 
its  best-fit  distance  of  d  ~  35  pc,  its  high  proper  motion  of 
1/75  yr-1  makes  it  a  member  of  the  halo,  with  an  estimated 
total  age  of  ^11  Gyr.  If  its  distance  is  d  <  13.5  pc,  it  could 
be  a  halo  or  a  thick  disk  object.  A  parallax  measurement 
for  J 1 102+41 13  is  needed  to  settle  the  question  of  its  origin 
and  to  determine  its  mass.  To  improve  the  data  available 
for  fitting  models  to  J 1 102+41 13,  deeper  u  photometry  and 
photometry  at  A  >  2.5  /x m  might  be  useful,  as  might  optical 
spectroscopy  with  a  higher  signal-to-noise  ratio  (S/N).  As  for 
the  models  themselves,  improvements  to  mixed  H-He  models 
are  known  to  be  needed.  In  addition,  better  H2  CIA  opacities 
(especially  at  optical  wavelengths)  and  more  examples  of  WDs 
at  reff  &  4000  K  are  needed  to  determine  how  well,  on 
average,  pure-H  models  at  such  temperatures  match  real  WD 
atmospheres. 

In  principle,  the  SDSS  can  be  used  to  place  an  upper  limit 
on  the  surface  density  of  red,  high-proper-motion  WDs  like 
J1 102+41 13.  In  practice,  additional  spectroscopy  beyond  that 
obtained  by  the  SDSS  will  be  needed  to  place  such  a  limit, 
because  objects  like  J 1 102+41 13  have  colors  too  similar  to 
those  of  the  stellar  locus  to  be  routinely  selected  for  SDSS 
spectroscopy.  (J  1 1 02+4 113  itself  barely  qualified  as  a  quasar 
candidate  based  on  preliminary  (“TARGET”)  SDSS  imaging 
reductions  and  no  longer  qualified  for  spectroscopy  based  on 
final  (“BEST”)  reductions.)  J1 102+41 13  is  the  only  g  <  19.7, 
g  —  r  >  1 .4  object  in  the  8417  deg2  SDSS  DR6  Legacy  imaging 
database  with  SDSS+USNO-B 1 .0  proper  motion  /x  ^  1 .5"  yr-1 
(Munn  et  al.  2004),  but  similarly  red  objects  with  /x  <  1"  yr-1 
that  lack  SDSS  spectroscopy  exist  in  the  database.  Whether  or 
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not  any  of  these  objects  turn  out  to  be  WDs,  J1 102+41 13  will 
remain  one  of  the  highest-proper-motion  cool  WDs  in  the  sky 
as  seen  from  Earth. 
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